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Abstract 



We analytically work out the long-term orbital perturbations induced by a homogeneous 
circular ring of radius and mass on the motion of a test particle in the cases (I): r > 
and (II): r < R^. In order to extend the validity of our analysis to the orbital configurations 
I ■ of, e.g., some proposed spacecraft-based mission for fundamental physics like LISA and 

ASTROD, of possible annuli around the supermassive black hole in Sgr A* coming from tidal 
disruptions of incoming gas clouds, and to the effect of artificial space debris belts around the 
Earth, we do not restrict ourselves to the case in which the ring and the orbit of the perturbed 
particle lie just in the same plane. From the corrections Aro^™"''''-' to the standard secular 
perihelion precessions, recently determined by a team of astronomers for some planets of the 
Solar System, we infer upper bounds on for various putative and known annular matter 
distributions of natural origin (close circumsolar ring with R^ = 0.02 — 0.13 au, dust ring 
with i?!- = 1 au, minor asteroids, Trans-Neptunian Objects). We find nii- < 1.4 x 10^** 
(circumsolar ring with R,- = 0.02 au), < 2.6 x 10~^ (circumsolar ring with R,^ = 0.13 
au), mr < 8.8 x 10~^ m® (ring with R^- = 1 au), rrir < 7.3 x 10~^^ Mq (astcroidal ring with 
Rr = 2.80 au), m,. < 1.1x10"" Mq (asteroidal ring with i?r = 3.14 au), < 2.0x10"^ Mq 
(TNOs ring with = 43 au). In principle, our analysis is valid both for baryonic and non- 
baryonic Dark Matter distributions. 

Keywords: Gravity - Asteroid Belt - Interplanetary Dust - Perturbation Methods 

1 Introduction 



The giant planets of the Solar System are surrounded by rings [Tj. They have widely been 
studied and have been the object of numerous scientific spacecraft-based missions [2l[3l[ll[5l[6l[7]. 
As we will see in more details in Section [21 there is the possibility that circumsolar massive 
rings also exist at small distances from the Sun [8l[9]. Moreover, the dynamical action of the 
ensemble of the minor asteroids of the asteroid belt between Mars and Jupiter can roughly be 
modeled as due to a continuous ring ^Q\. The same holds for the belt [HI [121 [E] of Trans 
Neptunian Objects (TNOs) [I4l[l5]. 

Recent observations suggest that similar structures exist in other stellar systems as well. 
[16] detected a ring-like structure at about 70 au from the Sun- type star HD 207129 [TT] with 
the IR observations collected by the instruments PACS [18] and SPIRE [19] of the Herschel 
space observatory [20]. For studies on Kuiper belt-like structures at about 30 — 50 au or beyond 
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around other nearby Solar-type main-sequence stars, see, e.g., [2T] and [22]. In particular, [21] 
detected 10 exo-Kuiper belts with masses ranging from to a minimu of 0.3 m(j to a maximum 
of 48 TTT-ij by means of far-IR observations from the APEX telescope [21] . 

On the other hand, rings are common features in different astrophysical objects as well. [25] 
suggested that supermassive black holes located at the center of galaxies may be shrouded by 
super-Oort clouds of comets and asteroids. [26] performed numerical solutions to the problem 
of black holes surrounded by uniformly rotating rings in axially symmetric, stationary space- 
times. [27j reviewed the properties of gravitating discs around black holes. Moreover, annular 
structures occur also in several ring galaxies (or R galaxies); they are objects with approximate 
elliptical rings and no luminous matter visible in their interiors \28\ 129]. Furthermore, in some 
circumstances, as the result of interactions between galaxies, a ring of gas and stars is formed 
and rotates over the poles of a galaxy, resulting in the polar-ring galaxies |30] . 

In addition to annular structures of natural origin, in the last decades analogous matter 
distributions of anthropogenic objects formed around the Earth [311[32]. Indeed, as first pointed 
out by [33], collisions among artificial satellites would produce orbiting fragments, each of which 
would increase the probability of further collisions, leading to the growth of a belt of debris 
around our planet in the Low Earth Orbit (LEO) regioro after just decades. [35] showed 
that the growth of such a debris population is primarily driven by high collision activities 
around 900 — 1,000 km altitude. Radar-based observational campaigns were implemented to 
characterize the orbital debris environment [36]. Much more long-lived debris equatorial rings, 
with lifetimes of the order of millennia, may form in the Geostationary Earth Orbit (GEO) 
regiorll as well [37l [38l [39] ; see, e.g. Figure 1.1 in [Ml P- 1]- For a recent overview, see [40] . 

Self-gravitating toroidal fluid configurations (without a central body) in Newtonian gravity 
were analytically studied by [HI SSI [S], [lUIlS], [l6], and [17]. In particular, HH iS] set up 
an approximation scheme of uniformly rotating, homogeneous, and axisymmetric rings which 
turned out to be an extremely good approximation for thin rings; see [l8]. For various physical 
properties of flat and toroidal ring models, see [l9] and [50], who exploited the disk schemes 
by [51] 5 [52] and [53]. For other approaches, see [5l] and [55]. [56] analytically considered 
the problem of a uniformly rotating, self-gravitating ring without a central body in Newtonian 
gravity. A general relativistic toroid around a black hole was studied by [57]. [58] and [59] 
considered three-dimensional models for the gravitational field of rings in the context of general 
relativity. Uniformly rotating, homogeneous and axisymmetric relativistic fluid bodies with a 
toroidal shape were investigated by [60]. [61] considered uniformly rotating homogeneous rings 
in post-Newtonian gravity. Uniformly rotating rings in general relativity were investigated also 
by [62]. [63] studied the black hole limit of rotating discs and rings. 

As far as the orbital perturbations induced by such mass configurations, [6H p. 58] and [65^ 
p. 195] treated the problem of the attraction of a continuous circular massive ring, and showed 
that its potential can be approximated by a complete elliptic integral of the first kind |66j . 
Actually, Gauss had already solved the problem and he used it to state his famous averaging 
theorem: replacing a perturbing body by an equivalent continuous mass spread over its orbit 
does not change the secular effects while it removes the periodic terms of the perturbation |67| . 
See also [68] • [69] analytically treated the equivalent problem of calculating the off-axis electric 
field of a ring of charge. Numerical computation of the acceleration due to a uniform ring 
can be found in |70j . The perturbing acceleration on an inner particle due to various external 
axisymmetric mass distributions (e.g. different kinds of thin rings, uniform and non-uniform 

lower bound of 76 rrifj was inferred by [21] for the star HD 95086 [23| . 
■^Satellites moving in it have typically altitudes h < 300 — 1, 500 km [341 p. 2] . 
''Geostationary spacecraft orbit at h = 35,800 km |34l p. 4]. 



2 



disk, torus) was worked out by [7T], [72], [73], [71] in the framework of the Pioneer anomaly 
|71] . |75] analytically worked out the force function of a slightly elliptical Gaussian ring, and 
treated its generalization to a nearly coplanar system of rings as well. For some approximate 
calculations of the orbital perturbations on the elliptical motion of a test particle moving inside 
annular distributions, see, e.g., [TB] and [77]. The nature of test particle motion in the presence 
of a rotating ring of self-gravitating matter around a Kerr black hole was investigated by [78j . 
The perturbations by a distribution of stars in the Galactic Center were investigated by [79] . 

In this paper, we analytically work out the long-term, i.e. averaged over one full orbital 
period, perturbations of a test particle moving about a central body and acted upon by a cir- 
cular massive ring. More specifically, after having reviewed various baryonic and non-baryonic 
scenarios involving annular matter distributions in Section [J] we consider the potential induced 
by the material annulus in the general case of an arbitrary ring-orbit configuration (Section [3|). 
Then, we analytically compute the long-term variations of all the osculating Keplerian orbital 
elements of a test particle due to a circular ring of homogenous linear density. As a by-product, 
we preliminarily put model- independent, dynamical constraints on the mass of several annular 
matter distributions from the orbital motions of some of the planets of the Solar System. In 
particular, we use the latest determinations of the supplementary secular precessions Am^'^'^^^'^ 
of the planetary perihelia [80]. In applying our results to the Solar System, we assume that 
the perturbed orbital motions occurs in the ring's plane itself. Table [1] shows that it can be 
considered a rather reasonable approximation for those planets of the Solar System for which 
[80] determined Aw'^"'"'^^ by using the mean Earth's equator at the epoch J2000.0 as reference 
{x, y} plane. We consider both the cases in which the planet moves outside the ring (Section 

Table 1: Inclination I, longitude of the ascending node /2 and argument of pericenter uj of the 
planets of the Solar System for which |80j determined the supplementary perihelion precessions 
Avu^^^^^^ with respect to the mean Earth's equator and equinox at J2000.0. / and H determine 



the orientation of the orbital plane in space, while uj yields the orientation of the orbit in its 
plane [ST] . 



HI), and inside it (Section [5|). As far as the hypothetical annuli are concerned, a complementary 
approach which may be followed consists of explicitly modeling them and solving for a dedicated 
parameter in fitting the models to the observations. 

We remark that our analysis is not necessarily limited to the major natural bodies of the 
Solar System (and of other stellar systems). Indeed, annular massive distributions like, e.g., 
the minor asteroid ensemble have an impact on several proposed space-based missions for fun- 
damental physics like ASTROD [82l [83l [84] , LISA [85l |M1 [87l [88] , GAIA [89l[90] relying upon 
spacecraft in Solar orbits of about 1 au. Moreover, man-made space debris belts orbiting the 
Earth may affect the orbital motion of terrestrial artificial satellites of interest for fundamental 
physics as well. Section [6] summarizes our findings. 

As a final remark, we point out that the present work has no pretensions of exhaustiveness 
about the whole ring problem. In particular, it purposely does not deal with the domain of 
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periodic orbits, where the potential in closed form expressions is used. For such a potential, 
there are orbits rather complicated, which likely could not be represented with only a few 
terms of the Legendre expansion used here. These orbits have a high variable orbital elements, 
and an averaging over the true anomaly like that performed in this paper gives no actual 
information. For studies on some general features of the dynamics around a massive annulus, 
like orbit stability for various configurations and the equilibrium of the system, see. e.g., 
[671 [m El [931 Ea Eg and [96] . 

Moreover, we do not pursue elegant and compact formulations which, if on the one hand 
may be of interest for the specialized reader, on the other hand may have little practical uses 
for the broad class of problems outlined here and for the allegedly broader audience interested 
in them. 

2 Baryonic and non-baryonic matter distributions in the Solar 
System 

2.1 Baryonic matter rings 

Is the Sun surrounded by a ring of rocky material lying at no more than a few Solar radii 
(1 Rq = 0.00465 au)? 

Such a hypothesis, along with that of a putative intra-Mercurial planet named Vulcan, was 
put forth for the first time by [97] in the form of a ring of asteroids, known also as "Vulcanoids" 
|98] , between Mercury and the Sun to explain the anomalous perihelion precession of Mercury of 
38 degrees of arc per century (arcsec cty~^) discovered by him in 1859. Later, [99], by including 
more observations in a new analysis of the motion of Mercury, came to an anomalous precession 

43 arsec cty ^. It was explained by [lOlj in terms of a zodiacal matter distribution whose 
gravitational attraction was equal to that of two ellipsoids, one inside Mercury's orbit and the 
other outside that of Venus. The hypothesis of a matter ring was supported also by [102j in 
a series of lectures given in 1906-1907 on the limits of the Newton's law of gravitation. As it 
is well known, [103j finally explained the anomalous perihelion precession of Mercury within 
the framework of its general theory of relativity which does not contain any ad-hoc, adjustable 
parameter. The possible existence of a circumsolar belt of small asteroids at about 0.1 au was 
revamped by [104| as a consequence of a reanalysis of some photographic plates taken during 
an eclipse in 1970. For more details, see |105j . [106| pointed out that Vulcanoids not smaller 
than 100 m may exist in a gravitationally stable band inside the orbit of Mercury, at distances 
of 0.06 — 0.21 au from the Sun, with an inclination to the ecliptic not exceeding 10 deg. See 
also [1071 11081 1109] . Their existence would be in agreement with the predictions of the scale 
relativity theory [llOj . [lllj . in their unsuccessful attempt to directly detect Vulcanoids in 
coronagraph images from the instrument LASCO ill2j carried onboard the spacecraft SoHO 
[113j . placed an upper bound on the size of such hypothetical objects which should not exceed 
60 km. In a previous SoHO/LASCO-based investigation, [109j concluded that the present-day 
population of Vulcanoids larger than 1 km should amount to less than about 1,800 — 42,000 
objects. 

The discovery of rings around some of the outer giant planets of the Solar System in the 
end of 70s [21 |3l H] enforced the possibility that a circumsolar massive ring exists as well [8] . 
Various evolutionary, physical, chemical and observational considerations have been provided 
about such a hypothesis. According to [8], the constituents would be > 10 km size boulders 

''its currently accepted value is 42.98 arcsec cty~^; see [100] for a discussion. 
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made of refractory material like graphite. As far as the radius i?r of a permanent ring structure 
around the Sun is concerned, the values [8] 

R, = 0.02 au = 4.3 Rq, 

(1) 



Rr = 0.13 au = 27.9 R, 



■05 



were proposed. The total mass would be much less than nij. ^ 6 x 10^^ kg |8]. [1141 1115] theo- 
retically predicted the existence of a circumsolar dust ring close to the Sun; for a review on the 
formation and possible observation of a dust ring, see |116| . Direct searches for a hypothetical 
circumsolar ring by exploiting its electromagnetic emission in the infrared (IR) portion of the 
spectrum were conducted yielding contradictory results |117j . especially as far as the nature 
and the dimensions of the putative constituents of the ring are concerned. Earlier observations 
of the F-corona brightness enhancement near 4 Rq by jll8j and [1191 1120] supported the ring 
hypothesis. |121] attempted to detect the existence of such a ring structure utilizing the IR 
observations taken during optimal viewing conditions of the total Solar eclipse of 16 February 
1980. The near-IR observations by |122] . who exploited the 1983 Solar eclipse, supported the 
existence of a circumsolar ring of dust lying approximately in the ecliptic plane at about 4 Rq. 
|123] . in their search for excess IR emission in the Solar equatorial plane during the eclipse of 
11 July 1991, concluded that the circumsolar dust ring may be a transient feature, perhaps 
due to the injection of dust into near-Solar space by a sun-grazing comet. jl24] argued against 
the existence of a circumsolar dust ring. According to [125] , the ring's constituents must be 
fine particles, instead of 10 km sized boulders at 4 Rq from the Sun. |126] concluded that it 
is not reasonable to suggest the presence of Si nanoparticles in the vicinity of the Sun. The 
hypothesis of a massive circumsolar ring made of macroscopic chunks of matter was recently 
reanalyzed by [9] who showed that, while a possible Solar ring structure at 4 Rq would not be 
stable over the age of the Solar System because of the Alfven drag |127] produced during even 
moderate Solar magnetic storms, a ring at ~ 27 Rq would, instead, be stable even during very 
intense Solar magnetic activity. 

|128] found evidence for a dust ring associated with the orbit of Venus in the data from 
the Helios-2 probe |129] . |130] and |131] pointed out the possible existence of a circumsolar 
ring of asteroidal dust with iZi- ~ 1 au; the azimuthal structure of this ring was predicted to 
be asymmetric, with the region trailing the Earth being substantially more dense than that in 
the leading direction. |132j . who analyzed data collected by the instrument Diffuse Infrared 
Background Experiment (DIRBE) [133] carried onboard the COsmic Background Explorer 
(COBE) satellite [134] . confirmed the existence of such a ring. [135] further investigated the 
structure of the Earth's circumsolar dust ring. For a review about dust bands in the Solar 
System and in extra-Solar planetary systems, see [136] . 

Other matter distributions, residing far from the Sun, which can be treated with the ap- 
proach outlined here are the ring of minor asteroids between Mars and Jupiter [10], and the 
TNOs' belt [m [121 (HJ- The action of the minor asteroids was modeled as due to a static 
circular ring at 2.8 au [TOl \m\ [IM [Ml EZ] or at 3.14 au [HOI [HH [HSl [l5]. [H [E] modeled 
the TNOs belt as a circular massive ring with R^- = 43 au. 

2.2 Non-baryonic Dark Matter distributions 

Until now we only considered rings made of ordinary baryonic matter. In principle, our analysis 
can extend to non-baryonic Dark Matter (DM) as well. 
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Its existence was postulated long ago to explain the discrepancy between the observed 
kinematics of some components of astrophysical systems like clusters of galaxies [143] and 
spiral galaxies [1441 1145j , and the predicted one on the basis of the Newtonian dynamics and 
the matter directly detected from the emitted electromagnetic radiation (visible stars and gas 
clouds). |146| postulated the existence of undetected, baryonic matter; today |147| it is believed 
that the hidden mass is constituted by non-baryonic, weakly interacting particles in order to 
cope with certain issues pertaining galaxy and galaxy clusters formation [148] . the Cosmic 
Microwave Background (CMB) [149j and primordial nucleosynthesis |150j . On cosmological 
scales, DM accounts for about 23% of the mass-energy content of the Universe [151j . A widely 
popular generic class of new-particle candidates is the Weakly Interacting Massive Particle 
(WIMP) scenario [1521 1153] . In it, DM annihilates with itself and interacts with the rest of the 
Standard Model (SM) via the weak interaction. The WIMP is typically defined as a stable, 
electrically neutral, massive particle which arises naturally in super symmetric SM extensions 
|154] . A pair of WIMPs can annihilate, producing ordinary particles and gamma rays. 

2.2.1 WIMPs 

DM particles in the Galactic halo, traversed by the Solar System during its Galactic journey, can 
end up on orbits gravitationally bound to the Solar System via purely gravitational three-body 
interactions with the planets |155j . However, when only gravitational interaction is considered, 
the inverse process of capture from the Galactic halo, leading to ejection of dark-matter particles 
from the Solar System, must be taken into account as well |155(ll56"tll57j . As a result, a detailed 
balance holds, so that a there is a rather small limiting density pdm of dark matter particles in 
the Solar System, not particularly larger than the mean Galactic value pbm ~ 4 x 10~^^ g cm~^ 
|158j . There are, however, different points of view on that issue according to which the situation 
with the true value for typical life-time of the captured DM in the Solar System would be far 
from being clear; see |159]ll60"lll61j . An increase in the population of DM particles bound to the 
Solar System can occur only if physical mechanisms other than pure gravitational scattering 
are at work like e.g., weak interactions with atoms in the Sun or the planets |162j . Within 
such a framework, |163j investigated a scenario in which WIMPs undergo weak scattering with 
atomic nuclei |164j in layers within the Sun close to its surface in such a way that they lose just 
enough energy to travel along closed Earth-crossing orbits with high eccentricities. Following 
such trajectories crossing the Solar surface, most of WIMPs are doomed to experience a second 
scattering event after just 10^ — 10^ revolutions which would significantly reduce their semimajor 
axes aw: they will end up in the Sun's core where they will ultimately annihilate with each 
other. However, if the layer of Solar matter traversed is quite thin, the lifetime of the grazing 
WIMPs gets larger. Moreover, during such an increased time the WIMPs experience small 
orbital perturbations by other planets of the Solar System which, through the |165j mechanism, 
stabilize them in stable orbital configurations at 1 au which do not intersect the Sun anymore 
over times comparable to the Solar System's lifetime |166j . More precisely, |166j obtained that 
the long-term survival of such WIMPs orbits is greater than 4.5 Gyr if their semimajor axes 
are smaller than half the semimajor axis of Jupiter, i.e. for aw < a7|_/2 = 2.6 au. |167j 
further studied such a scenario: anyway, according to [167| . the overall effect on the bound 
DM density should be small. The presence of such a stable WIMPs distribution at the Earth's 
location would have a great impact for the laboratory-based experimental searches of DM like, 
e.g., CDMSI [168], CDMSII [EH], DAMA/Nal [170] and its successor DAMA/LIBRA [TTT] . 
XENONIO [172] and ZEPLIN III [173]; for a recent introduction to dedicated DM experiments, 
see |174] . Thus, it is of great importance to put model-independent, dynamical constraints on 
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the total mass of such a putative distribution of DM at 1 au. 

2.2.2 Mirror matter 

Other possible stable matter distributions made of another kind of DM, the mirror matter |175t 
[T761 [T771 [T781 [T79l [THOl fTSTl [T82l [T83l [TM [T85l [T86] , have been hypothesized. It arises if instead 
of (or in addition to) assuming a symmetry between bosons and fermions, i.e. super symmetry, 
one assumes that nature is parity symmetric. For a popular overview of such a scenario, see 
[T87] . Mirror matter was supposed to exist in the Solar System; see [IMl IMl [HQl [Ml [T92l 
[1931 [TM] . 

2.2.3 CUDOs 

Recently, a new class of exotic objects in the Solar System has been proposed: the COmpact 
Ultra Dense Objects (CUDOs) |195| . They would be gravitationally self-bound objects made 
of non-baryonic (dark) matter fermions of mass 250 GeV to 100 TeV, or from massless fermions 
hidden by vacuum structure of similar energy scale |195| . The density at the center of a typical 
CUDO is more than 10^^ times higher than nuclear density, so that such extremely compact 
objects are small both in mass (sub-planetary, Mq < 0.1 M^)) and radius {Rc < 10 cm), making 
them difficult to observe by conventional astronomical methods |195j . Impact signatures in 
Solar System bodies have recently been proposed to constrain CUDOs' local abundances |196j . 
There may exist also CUDOs of other types, such as supersymmetric Q-balls |197[ I198[ I199j 
or CUDOs composed of standard model particles |200[ I201| . Q-ball collisions with Earth have 
been considered by j202j . 

3 Calculation of the potential of the ring for a generic orienta- 
tion of its plane in space and with respect to the orbit of the 
test particle 

For the sake of generality, in the following we will consider an arbitrary configuration of the 
ring-particle system. Indeed, in several astronomical and astrophysical scenarios involving, e.g., 
binary pulsar^ exoplanets, etc. the coordinate systenH adopted does not generally use any 
characteristic physical and/or orbital plane of the bodies of interest. Moreover, it is reasonable 
to expect that in several practical situations it cannot be assumed that the perturbing ring 
structures share the same plane of the perturbed bodies investigated. A striking example of 
non-coplanarity of several test particles orbiting their common primary is given by the stellar 
system orbiting the supermassive black hole in Sgr A* at the center of the Milky Way |205j : 
see [79] for a calculation of the orbital perturbations experienced by a S-type star due to a 
distribution of bodies. Recently, [206] discovered a dense gas cloud with rric = 3m® which is 
rapidly approaching the accretion zone of Sgr A* on a highly elliptic orbit (cc = 0.94) with an 
orbital period of 137 yr: it should reach the perinigricon in 2013. Over the past three years 
the cloud has begun to disrupt, likely by the tidal forces of the black hole [206] : it might, 
thus, create an annular matter distribution perturbing the S-type stars orbiting Sgr A*. Other 

^It is customarily adopted a system with the reference z axis oriented along the line-of-sight, from the object 
to the observer, and the reference {x, y} plane coinciding with the usual plane of the sky which is tangential 
to the celestial sphere at the position of the object. In it, the reference x axis is directed towards the North 
Celestial Pole (see, e.g., [2D31 p. 287] and p. 446]). 
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examples are given by the spacecrafts of the LISA [87] and ASTROD |82| constellations, whose 
orbital planes are planned to have quite different nodes Q (see, e.g., |207| ). In the case of the 
space debris belts surrounding the Earth [50], the GEO ring lies in the equatorial plane [57] . 
contrary to the orbits of several satellites potentially of interest for fundamental physics and/or 
other scopes which supposedly are perturbed by it. 

Let us choose a generic inertial frame endowed with a Cartesian coordinate system having 
the central body of mass M located at its origin O. Let r be the position vector of the perturbed 
test particle, while R denotes the position vector of a mass element draj- of the annulus. In 
general, the annulus is inclined by an angle /r to the reference {x, y} plane, and its intersection 
with the latter one is displaced by an angle with respect to the reference x axis. Thus, and 
]?r are the same for all the mass elements drrir of the annulus. Their position along it is picked 
up by an angle (/Jr, measured in the annulus's plane from its intersection with the reference 
plane. Thus, each dm^ has its own ip^. By assuming a homogeneous matter distribution, 
it will be drrir = prRdif^, where pr is the constant linear mass density. 

Quite generally, the Cartesian components of r can straightforwardly be expressed in the 
usual way [811 p. 51] in terms of the longitude of the ascending node the inclination / to 
the reference {x, y} plane, the argument of pericenter uj, the true anomaly /, which is the time- 
dependent, fast variable giving the instantaneous position of the test particle, the semimajor 
axis a, and the eccentricity e of the Keplerian ellipse of the test particle. 

Likewise, by assuming that the matter annulus has, in general, an elliptic shape, for the 
components of R we can write 

X = i? (cos /?r cos (/?r — COS /r sin /2r sin (/^j. ) , 

Y = i? (sin J?!- COS (/9r + COS Ir COS J7r sin , (2) 

Z = i?sin/i. sin^^r- 

Note that, in general, R can vary along the massive annulus according to a relation analogous 
to the standard expression of the Keplerian ellipse. Here we will assume R = throughout 
the ring. 

The potential exerted by a mass element dm^ in a point located at r is 

Gdnir Grrir dtpr 

dUr = — 1 =r- = :, (3) 

\r-R\ 27rg y/1 + a2 _ 2ar] {p,) ^' 

where 

(I) g = r, a = ^ for r > Rr, 

(4) 

(II) q = Rr, a = for r < R^, 

and 

r]{pr) = r-R (5) 
is the cosine of the angle between r and R. It explicitly depends on Lfj-. It can be written as 

T] (fr) = C cos ipr + S sinifj. , (6) 

with 

C = cos A/2cosii — cos/sin A/2sinn, 

(7) 

S = cos / sin A/2 cos n + (cos /cos /r cos A J7 + sin /sin /r) sin n. 
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with u = uj + f. From the point of view of the computation of the potential of the ring, 
I, n,u},f,Ir, fir entering eq. ^ are to be considered as fixed parameters. 
By expanding eq. ^ in powers of a < 1, it can be obtained 

dC/r = (I'd + ^'l + 1'2 + ^3 + ^4 + •••)c^V'r, (8) 

with 







2-Kq '■ 



^1 2-Kq 



*2 = ^{l-W)a\ (9) 



^3 = 



47rg 



^4 = 



The contribution of each term of eq. Q can be treated separately. First, the overall potential 
has to be computed by integrating over ip^ 



/•27r 

/ ^j{^r)d^r, j = 0,1,2,... (10) 
Jo 



We have 



1 



Uo = - 

Ui = 0, 

U2 = (11) 

U3 = 0, 



= -3GHk|l + :|(c4 + 54)-5[52+C2(l-|52)]} 



Then, the perturbations of the Keplerian orbital elements of the test particle induced by 
(t) , J = 0, 1, 2, • • • can straightforwardly be worked out with, e.g., the Lagrange planetary 

equations |8H pp. 251-252]. To this aim, the averages over one full orbital period must be 

worked out as ^ 

^^3) = (;^) / ' (*) J = 0, 1, 2, ... , (12) 

where Uj, j = 0,1,2,... are to be evaluated onto the unperturbed Keplerian ellipse of the 
moving test particle; as such, they are to be thought as functions of time t. Using the eccentric 
anomaly E \81\ p. 32] as fast variable of integration in eq. (I12p turns out to be computationally 
more convenient when q = since Uj oc , j = 1,2, . . .. 

Instead, when q = r and Uj oc r^^^^^\ j = 0,1,2,..., the use of the true anomaly / is 
computationally more suitable. 

Performing the integrations of eq. (|12p is, in general, quite cumbersome even for small 
values of j in eq. ([9]) . 
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If the annulus and the test particle's orbit are coplanar, the calculations are simpler. The 
potential of the ring can be worked out by using the standard mathematical result, valid for 
a < 1, 



2n 



27r Jq 1/1 + a2 - 2a cos ^ 



1 

— < 



K 



4a 
— Q 



K 



+ 



4a 
1 + Q 



where K{k) is the complete elliptic integral of the first kind [66]. Thus, the potential is 

" r , '^'"'f- . ^ ) (14) 

27rg 7o Vl + a2_2acos¥^r q \ i 64 J 

|208] yielded and expression in closed form, in terms of complete elliptic integrals. See also 
|209| I210| . and [211] who treated the equivalent problem of the computation of the averaged 
Hamiltonian in the case of a elliptical orbit perturbed by a circular perturber. It can be shown 
that eq. (jl4p yields just the acceleration at r obtained by |139j . 

4 The inner ring: r > Rj- 

In the case (I) of eq. (jH), it is 

{Uo) = (15) 
a 

Thus, provided that mr is by orders of magnitude less than the mass of the central body, all the 
Keplerian orbital elements are left unaffected, apart from the mean anomaly Ai. The Lagrange 
equation for its variation yields 

dM 2Gm, 

— — -nb = 3-. 16 



In it, nb = \/ GM / is the Keplerian mean motion connected to the orbital period by Tib = 
27r/Pb. 

The next non- vanishing averaged potential is {U2), which yields the following non-zero long- 
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term variations 



/ dl \ 3Gm,. / Hr \ 2 

\Tt/" 4nb(l-e2)2a3 \~ ) 

jO — sin Ji- (cos / cos Zi- + cosAf?sin7sin/i-)sinAf2, 

~dt)^~ 16nb (1 - e2)2 „3 \~) 

hO = — 2 cos 21 cos Afl CSC 7 sin 27, + cos 7 |^cos 27, (3 + cos 2Ar2) + 2 sin^ Afij , 
d-uj \ 3Gmr / 7?r 



dt / 64nb (1 - e2)2 V a 



PO, (19) 





'li 


'^7)+5sin| 






V2 , 






— go, 




a J 





pO = (3 - 4cos 7 + 5 cos 27) (1 + 3 cos 27,-) + 8 (3 + 5 cos 7) cos 2Ansin^ sin^ 7,-+ 

+ 2 cos An sec ^ — 

' dM \ 3Gmr , . 
-7r/~"b = ,,3/9 , go> (20) 



= + w is the longitude of pericenter. As expected, all the perturbations vanish in the 
limit iir — )• 0. If the ring and the planetary orbit share the same inclination to the reference 
{x,y} plane, but not the same orientation, i.e. for I = 1^-, fl ^ fl^, in general, non-vanishing 
perturbations depending on / and Q — Qt- occur. In the case fi = i7r, I 7^ Ir, the secular 
precession of the inclination / vanishes, while the rates of /?, ro, do not vanish, being functions 
of /, If. Finally, if both the ring and the orbit lie in the same plane, i.e. for 1 = 1^, fi = fl^ (cfr. 
Table then the inclination / and the node Q remain unaffected, while the other non-zero 
precessions are 



dw\ SGiTij- ^ 



dt / 4a3(l-e2)2nb 



(21) 



dM\ 3Gm, (R,V 



dt / 4a3(l-e2)3/2 



It can be noticed that eq. (j2ip and eq. (j22p are independent of the specific values of / and Q. 

In the case of (C/4), the resulting long-term rates of variation of the osculating Keplerian 
orbital elements for a generic ring-orbit configuration cannot be displayed in full since they are 
quite cumbersome; they are exact in the sense that no approximations in e, /, Jr, /?, /2r were 
adopted. In general, they do not vanish, and explicitly depend on / and A/2. For different 
inclinations but equal displacements with respect to the reference x direction (I ^ 1^, H = H^), 
it turns out that all the long-term rates of change are non-zero, explicitly depending on / and 
Jr. In particular, the shifts in e and / are of order 0{e) and O(e^), respectively. The case in 
which the inclinations are the same and the planes are shifted apart (/ = /r, i? 7^ H^) is more 
complicated. 
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Much simpler expressions are obtained for 1 = 1^, Q = Qt-. 

da\ , , 



dt I 



dl 
It 



(25) 



df2\ 

^>=0, (20) 



dw\ _ A5Gm, (l + fe^) /iJ^. 



/ 32a3(l_e2)4nb V « 



(27) 



^b = 77TTT^ ^777^1- ) • (28) 

nh 



4 



dt I 128a3 (1 - e2)^/2 

Also in this case, the non-vanishing rates neither depend on / nor on Q. 



4.1 The circumsolar ring with Ry. = 0.13 au 

For the sake of simplicity, we will assume coplanarity between the ring and the perturbed 
planetary orbit; eq. (|2T]) and eq. (|27|l yield 



\dt / 4a3 (1 - e2)2 



\ 15 (l + |e^) (rX 
8 (l-e2)2 V « / 



(29) 



In the case of Mercury (a^ = 0.387 au), and for Rj. = 0.13 au [9], eq. ([29|) yields 



= mr 3.08 X 10^2° mas cty~^ kg~\ (30) 

where mas cty~^ is a shorthand for milliarcseconds per century. An upper bound on nij. can 
be obtained from the latest determinations Aw!^'^^^ of the maximum allowed value for any 



unmodeled effect impacting the Hermean secular perihelion precession recently inferred by [80] 
from observations. [80], who did not model the Solar ring, obtained 

^^^meas) ^ ^ ^ ^ ^ g ^^^-l^ ^^^^ 

Actually, there is also a general relativistic effect, caused by the angular momentum of the Sun, 
which was not modeled by [80]: it is the gravitomagnetic precession of the perihelion |212| . 
which is predicted to be as large as 

Tjj^^'^^ = -2 mas cty"^ (32) 
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by using the value of the Solar angular momentum by helioseismology |213| . In the case of 
Mercury, also the impact of the mismodeling in the Sun's even zonal harmonic J2, modeled by 
[8Uj . must be taken into account. By assuming a 10% uncertainty in it |214| . it is 

w^^^^'^ = 2.5 mas cty"^ (33) 

Thus, both eq. (j32p and eq. (|33p must be subtracted from eq. ()3ip in inferring an upper bound 
on rrir. As a result, from 

0.5 mas cty"^ + rur 3.08 x 10"^° mas cty"^ kg"^ < 1 mas cty"^ (34) 

we have 

m, < 1.6 X 10^^ kg = 4.3 x lO"*' m®. (35) 

It is 3708 times smaller than the upper bound by [9]. 

One may legitimately wonder if the use of eq. (129p is accurate enough, given the relatively 
small distance of the Hermean orbit from the putative perturbing ring considered. Actually, 
the answer is positive since eq. ([35]) . inserted in eq. (p7|) . yields just 0.57 mas cty~"^. Thus, we 
conclude that the truncation level chosen in the expansion of the annulus' potential yielding eq. 
(|29p is fully adequate to deal with Mercury and a hypothetical perturbing massive ring with 
i?r = 0.13 au, in view of the present-day level of accuracy in constraining the supplementary 
Hermean perihelion precession. 

4.2 The circumsolar ring with Rj. = 0.02 au 

The same approach for Mercury and = 0.02 au yields 

m, < 8.3 X 10^° kg = 1.4 x 10"^ m©; 

it is tighter by about one order of magnitude than the upper limit proposed by [9] 
eq. (136]) and eq. (f27l) yield just 0.01 mas cty~^. 

4.3 The ring with = 1 au 

The perihelion precession of Mars {a^ = 1.524 au) can be used to effectively constrain the mass 
of a ring with i^r = 1 au. [80] obtained 

^^(meas) ^ ^ ^ -^^ ^^^-1 ^3^^ 

for the Arean supplementary perihelion precession. In the case of Mars, eq. (I29p and eq. (I37p 
yield 

rrir 2.1 X 10"^° mas cty"^ kg"^ < 0.11 mas cty~\ (38) 

which implies 

m, < 5.3 X 10^^ kg = 8.8 x 10"^ m®. (39) 

Since the Arean Lense-Thirring precession amounts to just '^^^^^ = —0.03 mas cty~^, i.e. 5 
times smaller than the present-day uncertainty [80] in determining the Arean supplementary 
perihelion, we did not include it in eq. (j38p . The same holds for the mismodeling in the Sun's 
quadrupole moment, since ro^"^^^ =0.01 mas cty~^. 

Also in this case, it turns out that the level of truncation chosen is adequate, given the 
current level of accuracy in constraining the supplementary Arean perihelion precession. Indeed, 
eq. ()39p . inserted in eq. ()27p . yields 0.05 mas cty~^, which is 3 times smaller than 0.15 mas 
cty-i. 



(36) 

. Note that 



13 



5 The outer ring: r < Rr 

In the case (II) of eq. the calculations are, in general, more involved. 

It turns out that, although {Uq) does not vanish, it does not contribute to the long-term 
variations of the osculating Keplerian orbital elements since it is just made of rrij- and -Rr- 
Indeed, it is Uq = —Gnir/Rr- Moreover, (f/i) = 0. As far as (C/2) is concerned, it contributes 
to the orbital changes with 



^ dt I 32»i|3 

El = cos 2ijj (-4 sin I sin 27, sin AB + 4 cos / sin^ /r sin 2Afi) + 

+ [(1 + 3 cos 27, ) sin + (3 + cos 27) cos 2Afi sin^ 7,-2 cos Afi sin 27 sin 27, j sin 2ij , 
\ 3Gm,- 



10 - 
12 - 



dn \ 

dt I 



(40) 



(210 + I2e^) . 



Srib^l - e2 7i3 

■ sin 7r (cos 7 cos 7, + cos A i2 sin 7 sin 7, ) sin A , 

■ (cos 7 cos 7r + cos Af^ sin 7 sin 7i-) [(3 + 5 cos 2ctj) sin 7, sin AT^ — 5 (cos 7, sin 7 — cos 7 cos AT^ sin 7i-) sin 2ci;] . 
3Gmr CSC 7 



(41) 



(200 - 



02c^ 



02 = 



(cos 7r sin 7 — cos 7 cos AJ^? sin 7,) (cos 7 cos 7, + cos Af2 sin 7 sin 7,) , 

(cos 7 cos 7r + cos Ar2sin7sin7r) [— ( — 3 + 5cos2ci;) (cos7rSin7 — cos 7cosAr^sin7r) — 5sin7i-sinAnsin2tt 
3Gmj 



(42) 



128 \/l - c^n^B^I 



(43) 



(PO + P2e^) , 

2 (1 + 3 cos 27,) ( 3 + 20cos2ijsin^ ^j)) + ''':°=2Afi (2-2 COS / + 5 cos 2a;) sin"^ Jr + 5 cos 2J ^1 + 3 cos 27r + 2 cos 2AJ7 cos luj sin^ /, 

"•(0) 

(^^^^ - 16cos2Afi (^5cos^ (^^^ co: 



40 sin 7 sin 27, sin Afi sin 2u + cos 7 4 (1 + 3 cos 27,) —1 + 5 cos 2ij s 



40 cos 2Afl sin j sin" 7, - 40 sin 7, sin 2Afi sin 2ii 



4 cos AfJ I 3 + 5 cos 27 + cos 7 (^6 - 20 cos" J cos 2uj ) | sin 27, tan 



P2 = -4 (1 + 3 cos 27,) -1 + 3 cos 7 + 10 cos 2u. 



OS 2uj + 3 sin I — I I sin 7,- + 
\2) I 



S COS Aft (-3 + 5 cos 2u) sec (^^^ sin 3 (^ sin 27, - 20 sec (^^^ (^3 sin ( ^ ^ ( ^ ) ) sin 27, sin A«sin2ij + 



- 80 cos ( — I sin 7, sin 2Aft sin 2aj , 
,2/ 



Gm, 



(mo + M2e^) , 



256nbi?i? 

= 4 (1 + 3 cos 27,) (7 + 21 cos 27 + 30 cos 2u sin^ 7) + 12 cos 2Afi |^10 (3 + cos 27) cos 2ij + 28 sin^ 7j sin^ 7,- 

- 48 cos Afi (-7 + 5 cos 2ij) sin 27 sin 27, + 480 sin 7 sin 27, sin Afi sin 2bj - 480 cos 7 sin^ 7, sin 2Aft sin 2ij, 

= 6 |2 (1 + 3 cos 27,) (1 + 3 cos 27 + 10 cos 2u sin^ 7) + 2 cos 2A« [lO (3 + cos 27) cos 2tj + 12 sin^ 7j sin^ 7, 

- 8 cos AJ7 ( — 3 + 5 cos 2aj) sin 27 sin 27, + 80 sin 7 sin 27, sin Af2 sin 2aj — 80 cos 7 sin" 7, sin 2Aft sin 2a;^ ; 



(44) 
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the semimajor axis is left unaffected. As expected, all the perturbations of eq. (jlO]) -eq. 
vanish in the limit -Rr — )■ oo. If the ring and the planetary orbit share the same inclination to 
the reference {x,y} plane, but not the same orientation, i.e. for 1 = 1,., Q ^ H^, in general, all 
the long-term changes of eq. (|40p -eq. (j44p do not vanish, depending on /, AJ7,a;; it is the case 
of, e.g., the three spacecrafts of LISA |2U7j . \i Q = Q^, / / Ij., then eq. (IlD|)-eq. dM]) are not 
zero, being functions of I,Ir,uj. Finally, if both the ring and the orbit lie in the same plane, i.e. 
for 1 = 1,., 17 = J?!- (cfr. Table [1]), then e, /, H remain unaffected, while the long-term variations 
of w and M. do not vanish, being equal to 



dw\ SVl — e^Gm,. 

(45) 



dt I 4nbi?? 
dM\ _ 7(l + fe^)Gm, 

^/-"^- ^ ^^^^ 

Note that eq. ([15]) and eq. ([16|) do not depend on the specific values of / and i?. 

The long-term variations of the Keplerian orbital elements due to (C/4) are quite involved 
for the general case I ^ I^, fi ^ fir- they cannot be explicitly displayed. Also in this case, no 
approximations in e,!,!^ were adopted. In the simpler coplanar case, i.e. for I = 1^, Q = H^, 
they are 

^\=0. (47) 



dt 

de 



(48) 



dQ 
It 



dw\ _ 45GmrVT^ (1 + |e2) / a ^ ^ 



dt I 32nbi?? V^: 



(50) 



(51) 



dt I 32nbi?i? V^r 

Also in this case, the precessions of eq. (|51|) and eq. (|52|) are independent of / and 12. 
For other analytical calculation, see, e.g., |77] . 



(52) 
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5.1 The ring with = 1 au 

As far as the perihelion precession in the coplanar case is concerned, eq. 



dw 
IF 



1 + ¥ ( 1 + TeM I ^ 



and eq. ([^T|) yield 
(53) 



By proceeding with eq. (|3ip as in the case of the inner ring, a tight constrain can be 
obtained by using the perihelion of Venus (o^ = 0.723 au), for which [80] obtained 



Aw, 



(meas) 



0.2 ± 1.5 mas cty"^ 



(54) 



(LT) 

By neglecting the Cytherean unmodeled Lense-Thirring and mismodeled J2 precessions ro^ 



-0.3 mas cty rij^^"^^^ 



0.2 mas cty ^, eq. ([53j) and eq. ([5^ yield 
m, < 2.1 X 10^9 kg = 3.4 x 10" " 



ma 



(55) 



It must be noticed that, for Rj- = 1 au, the Cytherean bound of eq. (j55|) is weaker than 
the Arean constraint of eq. (j39p . Anyway, the use of eq. (|53p to deal with Venus and the 
hypothetical ring at i?r = 1 au is adequate, given the current level of accuracy in constraining 
the supplementary Cytherean perihelion precession. Indeed, eq. ([5T]) . computed with eq. ([55]) . 
yields just 0.6 mas cty"^. 



5.2 The ring of minor asteroids between Mars and Jupiter 

The overall mass of the ring with which the dynamical action of the minor asteroids orbiting 
between Mars and Jupiter is usually modeled can effectively be constrained by using the perihe- 
lion of Mars and the result of eq. (j37p for its supplementary precession [80] . The impact of the 
individual largest perturbers of the asteroid belt on the motion of the inner planets, especially 
Mars and the Earth, was studied by, e.g., [2T5 | [2T6| [2T7 1 [2T8| [2T9 | [220] . As far as their impact 
on some proposed spacecraft-based missions for fundamental physics are concerned, see, e.g., 
[22n [2221122^ . 

By adopting [lOl [l371 UMl UMl [ZZ] = 2.80 au, eq. ([SSD, computed for Mars, yields 

TTir 7.5 X 10^^"*^ mas cty""*^ kg""*^ < 0.11 mas cty""*^, (56) 

from which it can be inferred 

rrir < 1.46 x 10^^ kg = 7.3 x 10"^^ Mq. (57) 

Also in this case, we did not consider the Arean Lense-Thirring effect in eq. (I56p because of its 
smallness. 

Our result of eq. (i57l) is smaller than the values obtained by [10] {m^ = (5±1) x 10"^° Mq), 
[T37] {rrir = (6±l)x 10-11 Mq), [l38] (r^r < 2.4x10"!° M©), [EH] (m^ = (3.4±1.5) x lO'^ M©), 
and [142] {m, = (6± 2) x lO"" Mq). In the case of [TiOl Ell [I42l [ll] R^ = 3.14 au, we obtain 
from Mars 

rur < 2.33 X 10^9 kg = 1.1 x 10"" Mq. (58) 

|142j and [141j obtained m,. = (l.OibO.3) x lO"!*^ Mq, which is one order of magnitude larger than 
eq. ([58]). [HQ] and [E] obtained m, = (3.35±0.35) x 10"!° Mq and m, = (8.7±3.5) x lO'^^ Mq, 
respectively. 
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The use of eq. ([8]) and eq. Q for r < i?;., yielding eq. (05]) and eq. (jST]) in the coplanar 
case, might be questioned because of the closeness of the Arean orbit to the asteroid ring. 
Actually, it is not the it can be a-posteriori inferred from eq. (j57p and eq. (j58|) . 

Indeed, it turns out that eq. ([5T]) . computed with both eq. ([57|) (i?r = 2.80 au) and eq. ([58]) 
[Rj- = 3.14 au), yields just 0.03 mas cty~^, which is 5 times smaller than 0.15 mas cty-^ Note 
that the same conclusion holds also if we adopt the upper limits of the figures for mr obtained 
by other authors. Indeed, for = 8 x 10"^^ Mq [1^ and = 2.80 au, eq. ^ yields 0.11 
mas cty~^. Moreover, mj- = 1.22 x 10~^^ Mq and Rj^ = 3.14 au jl5j . inserted in eq. (j5ip . give 
0.07 mas cty~^. Thus, we conclude that our expansion to q^, giving eq. (|53]l . is fully adequate 
to treat the case of Mars perturbed by the ring of minor asteroids within the present-day level 
of accuracy in constraining the supplementary precession of the Arean perihelion. 

5.3 The belt of the Trans-Neptunian Objects 

Concerning the mass of the TNO belt, modeled as a circular ring with |14] [T5] i?r = 43 au, it 
can be constrained with Saturn (oj^ = 9.582 au) by using |80j 

and eq. ()53p . Also in this case, there is no need to take into account the Kronian Lense- 
Thirring and 6J2 precessions since they are far smaller than 0.65 mas cty-i {w'^^^ = -1 X 

10~^ mas cty~^, roj^"^^^ = 2 x 10"^ mas cty~^). Thus, from 

mr 1.98 X 10"^^ mas cty"^ kg"^ < 0.8 mas cty"\ (60) 

it turns out 

ruj. < 4.03 X 10^2 kg = 2.0 x 10"^ M© = 6.6 x 10"^ m® = 42.7 mceres, (61) 
where we adopted |224j 

mCeres = 9.43 X 10^° kg (62) 

for the dwarf planet Ceres. 

For a comparison, according to [14], the minimum mass of the TNO belt would be as large 
as 110 masses of Ceres, while the maximum mass of the ring is expected to be 100 times the 
minimum mass, i.e. 1.1 x 10^ w,Ccrcs- In other tests by [14], the maximum mass surpasses 
the minimum mass by 25, 50, and 75 times. Thus, the Kronian bound of eq. (|6ip is much 
more stringent than that by [14] . jl5j inferred for the mass of the TNO ring the smaller value 
rur = (4.98 ± 0.14) x 10"^ Mq = (105 ± 3) mceres- The bound of eq. §1^ is also much tighter 
than those obtained some years ago by [76] from earlier corrections Arij^™^^^^ to the standard 
precessions of the perihelia of the inner planets |225] . 

We remark that, also in this case, the level of truncation considered is fully adequate, given 
the current level of accuracy in constraining the Kronian perihelion precession. Indeed, even 
by using the largest admissible value rrir = 5.12 x 10~^ Mq by |15j . it turns out that eq. ()5ip 
yields just 0.008 mas cty~^. 
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6 Summary and conclusions 



We considered a test particle orbiting a central body of mass M lying at the center of a 
circular massive ring of radius R^. We analytically worked out the long-term variations of the 
osculating Keplerian orbital elements a, e, /, i7, ca, of the particle induced by the first terms 
in the expansion of the gravitational potential of the annulus. We considered both the case 
in which the particle moves outside (r > R^) and inside (r < R^) it. We did not restrict to any 
specific mutual orientation of the ring and the particle's orbital plane which, in general, neither 
coincide (/ 7^ Ir, Q / i?r) nor lie in the reference {x, y} plane of the coordinate system adopted 
(/ 7^ 0,/i- 7^ 0). In general, the resulting perturbations affect all the particle's orbital elements, 
apart from the semimajor axis a. They depend in a quite intricate way on the inclinations /, /r 
and the nodes J7, J7j. of both the particle and the ring, apart from the case in which they are 
coplanar {I = 1^., Q = Q^). 

By assuming coplanarity, we applied our results to some hypothesized and existing annular 
distributions like circumsolar matter rings (Vulcanoids), a ring at 1 au, and the rings with 
which the dynamical actions of the minor asteroids in the asteroid belt and of the Trans- 
Neptunian Objects are usually modeled. By comparing our analytical results to the latest 
observational determinations of the supplementary precessions ^w^'^^'^'^ of the longitudes of 
the perihelia w of some planets of the Solar System, we tentatively inferred upper bounds 
on the masses mr of such massive annular distributions. As a complementary approach which 
could certainly be followed, a standard covariance-based analysis would consist of processing the 
same observational data records by explicitly modeling the annular mass distributions whose 
mass one wants to determine/constrain, and treat it as a solve- for parameter to be estimated 
in the fit of the models to the observations. In any case, it should not be considered as the 
only legitimate method which must mandatorily be followed to obtain meaningful results. In 
general, the issue of potentially occurring mutual cancelation of various competing effects can 
meaningfully be inspected analytically as well, provided that one is able to work out proper 
expressions adequate to the accuracy level of the data. In this respect, it is not true that 
the effects caused by the rings are the same as for, say, the Sun's oblateness: indeed, their 
analytical expressions are quite different. Moreover, in most cases, several competing effects 
were actually modeled, so that they can impact the supplementary perihelion precessions only 
with their mismodeling, which can well be taken into account in evaluating their "residual" 
effect. We stress that the constraints obtained by us cannot be considered as unrealistically 
tight; on the contrary, they might be regarded as relatively "generous" since, in inferring them, 
we assumed just that the entire range of variation of the supplementary perihelion precessions 
is due to the dynamical action of the rings. We also tested the validity of such an assumption by 
checking in each case if competing mismodeled/unmodeled known effects like the Sun's J2 and 
the Lense-Thirring effect were actually below the uncertainties released for the supplementary 
perihelion precessions; whenever it was not the case, like for Mercury, we took them into account 
as well. Claiming a-priori that there might still be room, in principle, for the action of any 
other unknown force conspiring in partly canceling allegedly larger ring effects makes little 
sense. After all, the effect of any sort of "Russell teacups" may well creep into the the solved- 
for rings' masses estimated in a covariance analysis, if they really existed in Nature and they 
were not modeled at all. Moreover, we successfully checked that the level of approximation 
used for our formulas is adequate, given the present-day level of accuracy in constraining the 
supplementary precessions of the planetary perihelia. On the other hand, apart from the fact 
that our constraints are compatible with other, previously obtained bounds, it is certainly not 
strange that they are smaller then them, given that we adopted the latest ephemerides. 
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Our results are summarized in Table [2j 

Table 2: Upper bounds on the mass of different circular matter rings of natural origin, whose 
radii are listed in au, inferred from the corrections Aw^'^'^^^^ to the standard perihelion 
precessions of some planets of the Solar System by [50]. We used eq. ([22]) and eq. ([35]) by 
assuming coplanarity between the rings and the planetary orbits. The upper bound for the 
TNOs {Rr = 43 au) is equivalent to 6.6 x 10 — 42.7 w-ccres' 



Rr (au) 




(max) /, N 

nir (kg) 


m^™'^^ (me, Mq) 


Planet adopted 


0.02 




8.3 X 10^0 


1.4 X lO"'' m® 


Mercury 


0.13 




1.6 X 10^9 


2.6 X 10"^ m© 


Mercury 


1 




5.3 X lO^s 


8.8 X lO"'^ m® 


Mars 


2.80 (minor 


asteroids) 


1.5 X 10^9 


7.3 X 10-12 


Mars 


3.14 (minor 


asteroids) 


2.3 X 10^9 


1.1 X 10-11 Mq 


Mars 


43 (TNOs) 




4.0 X 10^2 


2.0 X 10-8 Mq 


Saturn 



In principle, our results are not limited just to baryonic matter distributions, being valid 
also for non-baryonic rings made of Dark Matter. Concerning putative circumsolar massive 
rings, which may, in principle, be made of non-baryonic Dark Matter as well, their mass is 
no larger than about 1 x 10-^ — 3 x IQ-^ m©, while the mass of a ring at 1 au is ^ 9 x 
IQ-^ 771,0. Moreover, due their generality, our findings can be applied to different astronomical 
and astrophysical scenarios like Earth satellites, exoplanetary systems, stellar orbits around 
supermassive black holes, etc. 
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